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ABSTRACT 

The  strong  need  for  a  robust  and  real-tune  transmission  classification  can  be  found  both  in  civil  and  non-civil 
applications.  In  this  paper,  we  propose  a  new  technique  based  on  modem  and  low-complexity  object  detec¬ 
tion  approaches.  The  transmission  scanner  uses  a  combination  of  Haar-like  and  technical  features  to  detect 
and  classify  different,  co-occurring  narrowband  transmissions  within  wideband  signals.  The  transmission 
scanner  is  evaluated  with  recorded  real  world  data  in  order  to  fulfill  real  conditions.  The  evaluation  shows 
that  this  system  perfonns  very  well  with  99.5%  accuracy. 


1.0  INTRODUCTION 

To  know  what  is  on  air  has  always  been  a  strong  interest  of  civil  and  non-civil  institutions.  The  civil  area 
of  research  has  the  aim  to  detect  white  spaces  in  frequency  bands  in  order  to  establish  secondary  usage,  as 
it  is  done  in  Cognitive  Radio  (CR)  [1,2].  In  contrast  to  this,  the  non-civil  institutions  try  to  identify  every 
signal  or  a  special  type  of  signal.  This  signal  interception  is  called  Communications  Intelligence  (COMINT). 
Both  benefit  from  increasing  instantaneous  receiver  bandwidths,  made  possible  by  Software-Defined  Radio 
(SDR)  [3].  The  wider  the  bandwidth,  the  higher  the  knowledge  gain  can  be.  The  drawback  of  such  a  wide¬ 
band  receiver  is  that  the  requirements  for  a  real-time  signal  processing  in  CR  and  COMINT  applications  are 
higher  than  before,  as  the  search  space  increases. 

In  Fig.  1  the  traditional  procedural  method  is  shown.  A  Wideband  Receiver  (WBR)  digitizes  the  antenna 
signal  with  a  specific  sampling  frequency.  The  typical  sampling  frequencies  for  the  High  Frequency  (HF), 
Very  and  Ultra  High  Frequency  (VUHF)  bands  vary  between  fs  =  625  kHz  and  fs  =  20  MHz.  The  dis¬ 
cretized  signal  is  led  to  a  Wideband  Signal  Detection  (WDET).  The  most  approaches  for  this  stage  are  based 
on  an  energy  detection.  They  use  no  a  knowledge  to  extract  the  information  regarding  the  center  frequency, 
bandwidth,  transmission  start  and  duration.  The  information  of  WDET  is  used  within  the  Digital  Down  Con¬ 
version  (DDC)  to  cut  the  narrowband  signals  out  of  the  wideband  signal  and  transfer  them  to  the  Narrowband 
Transmission  Classification  (NTCL).  The  NTCL  extracts  specific  features,  e.  g.,  modulation  type  and  baud 
rate,  in  order  the  get  the  information  about  the  used  transmission  standard,  e.  g.,  STANAG-4285  or  PACTOR 
II  [4].  Approaches  for  the  modulation  classifier  are  presented  in  [5-8].  The  evaluation  of  most  of  these 
algorithms  has  been  done  without  consideration  of  inaccuracies  of  the  first  stage  and  furthermore  almost 
all  test  signals  have  been  created  without  consideration  of  transmission  standards.  Therefore,  two  problems 
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ORGANIZATION 


Figure  1:  Processing  chain  with  separated  wideband  signal  detection  and  narrowband  transmission  classification. 


arise  in  real  applications.  First,  most  of  the  classification  approaches  used  in  NTCL  extract  features  or  tech¬ 
nical  parameters  which  loose  their  quality  in  case  of  detection  inaccuracies,  e.  g.,  frequency  offset.  Second, 
several  transmission  standards,  e.  g.,  STANAG-4285  [4],  change  the  modulation  index  or  the  type  during  the 
transmission,  or  use  a  double  modulation,  e.  g.,  PACKET-RADIO  [4].  Therefore,  a  feature,  which  is  able 
to  distinguish  specific  modulation  types,  might  fail  in  applications  where  a  signal  with  different  modula¬ 
tion  types  commonly  occurs.  The  biggest  problem  of  the  block  diagram  shown  in  Fig.  1  is  the  separtion  of 
wideband  and  narrowband  processing.  The  resulting  quality  of  signal  detectors  degrades  in  case  of  lower 
signal-to-noise  ratios  (SNRs)  or  stronger  fading  effects.  This  leads  to  incorrect  information  for  the  DDC 
block  and  to  useless  narrowband  signals  for  the  NTCL  stage.  In  most  of  the  cases  the  classification  results 
are  incorrect.  Thus,  it  is  advisable  to  avoid  a  separation  between  wideband  and  narrowband  processing  or  to 
combine  signal  detection  and  classification. 

With  tremendously  increased  processing  power,  more  research  has  been  done  in  SDR  in  the  last  years. 
Especially  in  the  area  of  SDR  and  CR  new  strategies  of  an  improved  combined  detection  and  classification 
emerged.  Beside  the  typical  energy  based  detectors  [9]  the  first  detectors  using  a  priori  knowledge  came  up. 
These  detectors  are  based  on  features,  for  example,  second-order  cyclostationarity  [10,  11]  or  higher-order 
statistics.  They  are  reliable  even  at  lower  SNRs,  but  the  computational  effort  is  very  high.  Within  the  pre¬ 
sented  literature  it  is  mostly  assumed  that  the  center  frequencies  of  the  signals  of  interest  are  well  known, 
which  is  true  for  most  so-called  cooperative  applications.  In  non-civil  or  non-cooperative  applications  the 
only  assumption  that  can  be  made  is  that  there  might  be  a  signal  of  interest  at  an  unknown  position.  This 
requires  a  scanning  classification  system  which  is  possible  to  perform  the  scanning  and  classification  process 
in  real-time. 

In  this  paper  we  present  a  flexible  combined  transmission  detection  and  classification  for  both  coopera¬ 
tive  and  non-cooperative  applications.  It  will  be  trained  for  a  common  HF  transmission  standard.  The  system 
is  evaluated  with  real  signals  in  order  to  show  the  robustness. 

The  paper  is  organized  as  follows:  Section  2.0  describes  the  principle  of  transmission  scanning.  In 
Section  2.1,  we  present  the  features  used  for  the  classification.  Section  2.2  analyzes  the  classification  system. 
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Section  3.0  reports  on  the  performance  of  our  system.  In  Section  4.0,  we  present  conclusions  along  with 
recommendations  for  future  improvements. 


2.0  TRANSMISSION  SCANNER  FRAMEWORK 

The  basis  of  our  approach  is  the  spectrogram  of  a  wideband  signal.  In  Fig.  2  only  a  small  part  of  a  wideband 
spectrogram  is  shown.  Many  experts  are  able  to  detect  and  classify  transmissions  only  by  looking  at  this 
special  visualization.  Thus,  the  principle  of  our  proposed  transmission  scanner  is  comparable  to  well-known 
object  detection  techniques  (e.  g.,  face  detection  [12]),  which  are  applied  to  images.  In  short,  the  scan  is  a 
series  of  detections  and  classifications  of  different  masks  which  are  moved  over  the  image.  In  our  case,  the 
image  is  comparable  with  the  complex  spectrogram  X 
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where  |x*|  are  the  bins  of  the  power  spectrum  and  ts  is  the  number  of  short-time  spectrums  used  for  classi¬ 
fication.  Xj  is  computed  within  the  wideband  receiver,  a  ComCat™-Tuner  delivered  by  MEDAV  [13].  With 
an  applied  FFT  length  /fft  =  4096  samples,  a  sampling  frequency  fg  =  625  kHz  and  a  Hann  window  with 
50%  overlapping,  we  receive  a  frequency  resolution  A /  =  153  Hz  and  a  time  resolution  A t  =  0.003  sec1. 
The  mentioned  mask  is  fitted  to  the  transmission  we  are  looking  for.  In  Section  3.0  the  transmission  of 
interest  is  a  STANAG-4285  signal.  It  is  a  continuous,  fixed  frequency  and  3.0  kHz  wide  signal.  It  became 
evident  that  0.500  s  signal  length  is  sufficient  for  a  robust  classification.  Thus,  the  mask  is  3.3  kHz  wide  and 
0.500  sec  long.  Fower  time  lengths  are  possible,  but  will  result  in  a  lower  classification  rate.  This  leads  to  a 
mask  Xm,  /m  =  22  bins  wide  and  =  168  bins  long, 


XM  (u,v) 


X  [ u ,  n] 


X[u  +  tM,w] 


X[u,v  +  /M] 

X  [u  +  tM,  v  +  /M] 
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where  u  =  [1, ...,  —  %]  is  the  position  in  time  direction  and  v  =  [1, ...,  /fft  —  /m]  the  position  in  fre¬ 

quency  direction  of  X.  In  Fig.  3  the  block  diagram  for  the  new  system  is  shown.  By  sweeping  of  Xm  or 
rather  varying  u  and  v,  all  possible  mask  positions  within  X  are  sent  to  the  wideband  feature  extraction 
and  then  to  the  wideband  classification.  The  wideband  feature  extraction  is  described  in  Section  2.1  and  the 
wideband  classification  in  Section  2.2. 

'These  resolutions  are  adapted  to  the  High  Frequency  (HF)  range.  For  the  Very  and  Ultra  High  Frequency  (VUHF)  bands  other 
resolutions  should  be  chosen. 
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Spectrogram  of  multiple  transmissions. 


Frequency  in  [kHz], 

(a) 


Spectrogram  of  multiple  transmissions. 


Frequency  in  [kHz]. 

(b) 


Figure  2:  Part  of  an  absolute,  logarithmic  wideband  spectrogram.  It  contains  several  transmissions,  for  example,  STANAG- 
4285  and  BAUDOT.  A  STANAG-4285  mask  XM  with  /M  =  22  bins  and  tM  =  168  bins  is  exemplarily  swept  over  the  spectrogram. 
In  Fig.  2(a)  XM  lies  not  on  a  STANAG-4285  transmission.  In  Fig.  2(b)  XM  lies  on  a  STANAG-4285  transmission. 
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Combined  Wideband  Processing 


STANAG-4285 


PACTOR II 


Figure  3:  Processing  chain  with  a  combined  wideband  signal  detection  and  transmission  classification. 


2.1  Wideband  Features 

For  transmission  classification,  it  is  desirable  to  define  iVpeat  different  feature  types.  There  are  features  which 
are  very  efficient  but  not  definite,  and  there  are  features  which  allow  exact  decisions  but  are  very  inefficient. 
Thus,  the  proposed  framework  gives  us  the  possibility  to  include  different  feature  groups  and  to  benefit  from 
the  different  advantages.  At  the  end  of  the  feature  extraction,  we  obtain  a  feature  vector  f 


f  =  [fl,  .  .  .  ,fjVpeJ 

with  N 


(3) 


N  =  Ni  +  ...  +  Nn Feat 


(4) 


dimensions. 

In  the  research  area  of  object  detection  optical  features,  like  Haar-like  features,  are  state  of  the  art. 
They  are  very  efficient  and  obtain  a  good  quality.  In  our  approach,  we  apply  the  Haar-like  features  to  the 
logarithmic,  absolute  mask  XM.iog 


XM.iog  =  10  •  log10  (|Xm|2)  ■  (5) 

The  original,  complex  mask  Xm  gives  us  the  possibility  to  extract  more  information,  which  is  invisible 
in  the  logarithmic,  absolute  spectrogram.  In  addition  to  the  optical  Haar-like  features,  we  use  technical, 
phase-related  features,  too. 
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2.1.1  Haar-like  Features 

The  evaluation  of  every  single  bin  within  the  mask  does  not  give  any  information  other  than  signal  power 
and  phase  information,  which  is  not  sufficient  for  the  classification  process.  So,  first  of  all  we  need  features 
that  encode  the  relation  of  different  signal  power  values  within  one  observed  mask  XM.iog-  Haar-like  features 
are  typically  used  to  describe  contrasts  within  images.  Applied  to  our  problem,  these  features  characterize 
the  two-dimensional  energy  distribution. 

Our  Haar-like  feature,  which  is  shown  in  Fig.  4,  set  was  inspired  by  Papageorgiou  et  al.  in  [14].  We 
reduced  the  number  of  Haar-like  features  to  a  radio  signal  specific  set.  It  contains  now  20  2-D  windows,  and 
can  be  used  for  different  transmission  types.  Depending  on  the  Haar-feature  that  is  used  for  extraction,  dif¬ 
ferent  tasks  can  be  fulfilled.  Fig.  4(f)  and  Fig.  4(g)  are  for  the  symbol  alternations  of  a  2-FSK  transmission, 
whereas  Fig.  4(h)  is  useful  for  the  bandwidth  estimation. 


(0 


(c)  (d)  (e) 


Figure  4:  Haar-like  features  used  for  detection  and  classification.  Dark-colored  areas  have  negative  and  light-colored  areas 
positive  weights. 


The  Haar-like  features  are  applied  to  different  positions  with  variable  scale  factors  within  the  mask  Xm, 
as  it  shown  in  Fig.  5.  In  this  case,  the  resulting  feature  value  is  the  difference  between  the  rectangle  A  and 
rectangle  B.  Considering  all  possible  positions  and  scaling  factors,  Nu aar  =  362.720  Haar-like  features  are 
extracted.  Viola  et  al.  showed  in  [12]  a  very  fast  computation  scheme  based  on  integral  images,  which  has 
been  adapted  in  our  implementation. 


2.1.2  Phase-related  Feature 

In  this  paper  we  introduce  a  technical  feature  which  is  especially  suited  to  detect  periodic  sequences  within 
a  transmission.  These  sequences  are  part  of  the  standard  and  often  used,  in  order  to  enable  equalization 
algorithms.  This  feature  is  based  on  a  cepstral  analysis,  which  is  suitable  for  the  detection  of  periodic 
values,  as  described  by  Bogert  et  al.  in  [15].  In  contrast  to  the  cyclostationary  feature,  presented  by  Gardner 
et  al.  in  [10],  the  cepstral  feature  is  modulation  type  independent.  Thus,  it  is  suitable  for  transmissions 
with  changing  modulation  types  or  double  modulated  content.  Additionally,  the  distance  between  training 
sequences  is  often  unique  for  different  transmission  standards.  The  cepstrum  qceps  (i)  of  a  complex  signal 
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A  Haar-like  feature  applied  to  the  mask  Xj\/. 
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Figure  5:  The  Haar-like  feature,  which  is  shown  in  Fig.  4(a),  is  applied  to  the  mask  XM  at  a  certain  position  and  with  specific 
scaling  factors. 


q  (t)  can  be  expressed  as  a  function,  as  follows 


QCeps  (t) 


log 


Q(w) 


~jwtdw 


(6) 


where  Q  (w)  is 


/+oo 

q  ( t )  e~jwtdt.  (7) 

-OO 

In  our  case,  q  (i)  is  the  downsampled  narrowband  signal,  dependent  on  the  position  of  the  mask  Xm  within 
X.  In  Fig.  6  the  cepstrum  of  two  transmissions  is  plotted.  It  shows  the  position  of  a  training  sequence  and  its 
harmonics.  Exploiting  the  harmonic  structure  within  the  cepstrum,  the  distance  between  the  sequences  can 
be  extracted  and  used  for  classification.  In  both  cases  a  peak  analysis  gives  the  correct  result  of  106.6  msec 
between  the  sequences.  In  Fig.  6(b)  it  emerged  that  the  cepstrum  is  very  robust  against  fading  and  low  SNR. 
With  iVcepst  =  1,  we  obtain  a  feature  vector  f  with  N  =  362.721  different  elements. 

2.2  Wideband  Classifier 

Given  a  set  of  labeled  feature  vectors,  which  are  computed  from  the  signal  and  rejection  patterns,  any  ma¬ 
chine  learning  approach  can  be  used  to  learn  a  functional  relationship  between  feature  values  and  classes. 
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Although  the  Haar-like  features  can  be  computed  very  efficiently,  extracting  the  complete  set  is  very  expen¬ 
sive.  So,  we  decided  to  use  the  Adaptive  Boosting  (AdaBoost)  algorithm,  which  is  well-known  in  object 
detection  and  first  published  by  Freund  et  al.  in  [16].  AdaBoost  trains  T  weak  classifiers  hi  (Xm)  repeatedly 
with  t  =  1 , ,T.  The  final  strong  classifier  is  a  weighted  combination  of  the  iterated  hypotheses  ht.  Gen¬ 
erally,  for  every  iteration  a  weak  classifier  which  evaluates  a  threshold  Qt  on  a  single  feature  ft  is  trained  to 
distinguish  between  the  signal  class  c\  and  the  rejection  class  c 2: 

h,  (XM)  =  |CI’  /‘(Xm,<6‘  (8) 

I  C2 ,  otherwise 

After  each  iteration,  the  training  set  is  reweighted  with  respect  to  the  error  rate  of  the  last  weak  classifier. 
As  T  is  chosen  to  be  smaller  than  the  overall  number  of  available  features  N,  the  number  of  features  will  be 
reduced  during  the  training  process.  N'  will  remain  after  the  training: 

N'  «  N  with  N'  =  A^aar  +  iV'epst.  (9) 

It  can  happen  that  Ar^aar  or  iV£epst  is  zero.  In  this  case  the  corresponding  feature  group  will  not  be  used  for 
the  classification  process.  For  more  details  regarding  AdaBoost,  see  [16]. 


3.0  EXPERIMENTAL  RESULTS 

The  task  is  to  detect  STANAG-4285  transmissions  in  wideband  signals  without  prior  knowledge,  e.  g.,  the 
exact  position,  and  reject  noise  and  all  other  transmissions.  For  the  experiments,  we  trained  three  different, 
monolithic  AdaBoost  classifiers  with  T  =  10  weak  classifiers.  The  first  classifier  is  trained  with  both  feature 
groups  fjjaar  and  fceps,  while  the  others  are  trained  each  with  feature  group  fnaar  or  fceps- 

Unfortunately,  there  exists  no  realistic  labeled  data  from  real-world  transmissions.  Thus,  the  wideband 
signals  used  for  this  work  were  acquired  and  manually  labeled  by  MEDAV.  Researchers,  interested  in  a 
comparison  with  their  own  approach,  can  have  access  to  the  data  upon  request.  The  wideband  signal  for 
training  has  a  sampling  frequency  fs  =  625  kHz  and  a  center  frequency  fc  =  8.350  MHz.  The  wideband 
signal  for  the  evaluation  has  the  same  sampling  frequency  but  a  center  frequency  of  fc  =  8.800  MHz. 
To  provide  typical  rejection  samples  for  this  scenario,  the  rejection  class  contains  examples  of  noise  and 
examples  of  common  HF  transmissions,  e.  g.,  SITOR-ARQ,  BAUDOT,  GW-PACTOR,  GW-OFDM,  J3E- 
USB  and  A1E  (described  by  Prosch  in  [4]).  The  used  wideband  signals  are  an  important  part  of  the  realistic 
evaluation,  as  they  represent  a  typical  HF  scenario.  They  contain  a  lot  of  unknown  transmissions  and  time- 
variant  wideband  noise.  In  Table  1  the  class  distribution  is  shown. 

Table  2  shows  the  results  of  the  evaluation.  Compared  to  the  trainings  with  only  one  feature  group,  we 
obtained  the  best  classification  rate  of  99.5%  for  real  world  scenarios  with  the  combination  of  fnaar  and  fceps- 
In  Fig.  7,  only  a  part  of  the  results  is  plotted. 

The  total  positive  rate  for  the  rejection  class  is  about  99.7%,  and  the  total  positive  rate  for  the  signal  class 
is  96.3%.  This  indicates  the  fact  that  the  combination  of  both  groups  is  necessary  for  a  robust  classification. 
A  high  accuracy  for  the  rejection  class  is  highly  appreciated  in  COMINT  applications,  because  the  post  pro¬ 
cessing  of  missclassified  signals  is  very  expensive. 
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Class 

Name 

Training 

Test 

Cl 

STANAG-4285 

4887 

1581 

C2 

NOISE 

12434 

4646 

C2 

SITOR-ARQ 

2898 

2844 

C2 

GW-PACTOR 

13386 

6960 

C2 

BAUDOT 

5544 

6336 

C2 

GW-OFDM 

9504 

5244 

C2 

J3E-USB 

2088 

7032 

C2 

A1E 

8730 

0 

Table  1 :  Number  of  patterns  for  training  and  test  set.  Due  to  the  realistic  properties  of  the  wideband  signals,  it  is  possible  that 
some  classes  have  more  or  less  samples,  e.  g.,  A1E.  Both  training  and  test  set  contain  transmissions  with  a  SNR  between 
7  dB  and  30  dB. 


Training  and  evaluation  for  T  =  10 

Feature  vector  fHaar,  fops 

TP  Rate 

FP  Rate 

Precision 

Recall 

F-Measure 

ci  =  STANAG-4285 

0.963 

0.003 

0.932 

0.963 

0.947 

c2  =  REJECTION 

0.997 

0.037 

0.998 

0.997 

0.997 

Weighted  Avg. 

0.995 

0.035 

0.995 

0.995 

0.995 

Feature  vector  [fHaar] 

ci  =  STANAG-4285 

0.822 

0.019 

0.671 

0.822 

0.738 

c2  =  REJECTION 

0.981 

0.178 

0.991 

0.981 

0.986 

Weighted  Avg. 

0.973 

0.171 

0.977 

0.973 

0.975 

Feature  vector  fceps 

ci  =  STANAG-4285 

0.917 

0.015 

0.745 

0.917 

0.822 

c2  =  REJECTION 

0.985 

0.083 

0.996 

0.985 

0.990 

Weighted  Avg. 

0.982 

0.080 

0.985 

0.982 

0.983 

Table  2:  Evaluation  of  the  classifiers  with  real  world  signals.  All  possible  combinations  of  the  feature  groups  fHaar  and  fCeps 
have  been  considered. 


Using  the  integral  feature  computation  and  some  other  optimizations,  we  could  greatly  reduce  the  com¬ 
plexity.  We  have  been  able  to  run  about  13.108  classifications  per  second2,  which  is  a  great  step  to  real-time 
CR  and  COMINT  applications. 


4.0  SUMMARY  AND  OUTLOOK 

The  paper  introduced  a  novel  method  for  transmission  scanning  based  on  well-established  object  detection 
techniques.  In  order  to  improve  the  accuracy  we  used  beside  the  typical  Haar-like  features,  also  a  new  cep- 
stral  feature.  A  combination  of  both  feature  groups  exploited  the  particular  benefits.  An  adaptation  to  new 
transmissions  or  to  VUHF  scenarios  is  possible  at  any  time.  With  signal  recordings  provided  by  MEDAV, 
we  could  show  that  the  system  performs  successfully  and  very  fast  under  real  circumstances.  We  obtained 
very  good  results  for  a  common,  real  HF  scenario. 

2  A  single  threaded  implementation  ran  on  a  Intel®Core™2  Duo  CPU  T8300  at  2.40  GHz. 
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In  future  we  will  implement  an  efficient  multiclass  version  of  the  transmission  scanner,  based  on  feature 

sharing.  In  order  to  show  the  flexibility  of  our  system,  we  will  focus  on  VUHF  transmission  standards. 
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Spectrogram. 
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Figure  6:  Spectrograms  and  cepstra  of  both  a  high  SNR  and  a  low  SNR  STANAG-4285  transmission.  The  dimensions  of  the 
spectrogram  are  equivalent  to  the  dimensions  of  XM. 
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Results  of  the  transmission  scanner. 


Figure  7:  Part  of  an  absolute,  logarithmic  wideband  spectrogram.  It  contains  the  results  of  the  transmission  scanner,  marked 
as  green  boxes.  This  classifier  has  been  trained  with  a  combination  of  fHaar  and  fceps-  Every  STANAG-4285  transmission  has 
been  found.  There  is  no  incorrect  hit  in  this  part  of  the  wideband  spectrogram. 
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